Partitioning optical solitons for generating entangled light beams 
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It is shown that bipartition of optical solitons can be used to generate entangled light beams. 
The achievable amount of entanglement can be substantially larger for A'^-bound solitons (A*' = 2, 3) 
than for the fundamental soliton (A*' = 1) . An analysis of the mode structure of the entangled beams 
shows that just A*' modes are essentially entangled. In particular, partitioning of the fundamental 
soliton effectively produces 2- mode squeezed light. 
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Quantum coherent dynamics of atomic and photonic 
systems has been a subject of increasing interest stimu- 
lated by the rapidly developing field of quantum informa- 
tion processing and quantum computing [Q-^. Intrinsic 
parallelism of quantum state evolution renders it possible 
to overcome the problem of exponential time required by 
a classical computer for solving complex problems such 
as integer factorization 1^,0, discrete logarithm and 
search H. High sensitivity of quantum coherent dynam- 
ics to eavesdropping can be used for developing secure 
encryption protocols realized in fiber optical systems 
& 

One of the most remarkable features of quantum coher- 
ence is entanglement, on which, e.g., quantum teleporta- 
tion is based . Entanglement can be regarded as being 
the nonclassical contribution to the overall correlation 
between two parts of a system (see, e.g., |l^Jl^). Typi- 
cally, discrete systems that are built up by qubits have 
been considered. Concepts that use continuous quan- 
tum variables have recently attracted increasing publicity 
[p^-p^ . An illustrative example is the teleportation of co- 
herent states by means of entangled light of the squeezed- 
vacuum type . 

A way to generate entangled light is to use two fields 
each of which is prepared in some nonclassical state 
and combine them at a linear four-port device like a 
beam splitter or a coupler. In particular, the entangled 
squeezed light used in is produced by combining 
two light beams independently squeezed via x^^^ para- 
metric processes. The possibility of realizing entangled 
pulse sources by combining two separately squeezed op- 
tical solitons are discussed in |Q . Alternatively, a nonlin- 
ear coupling between two fields can prepare them in an 
entangled (continuous-variable) state. For example, the 
T^(2) process of parametric down-conversion has been ex- 
tensively used as a source of entangled beams and pulses 
(see, e.g., ^ , p2| ). In this letter we show that bipartition 
of (fundamental and A^-bound) optical solitons formed 
in a medium with a Kerr- type x^'^-' nonlinearity yields 
entangled light fields, because of the internal quantum 
correlations of the solitons. 

From classical optics it is well known that the Kerr 
nonlinearity can compensate for dispersion-assisted pulse 



spreading or diffraction-assisted beam broadening (see, 
e.g., [^,g4j). In both cases, the evolution of the (com- 
plex) field amplitude a{x, t) can be described by the non- 
linear Schrodinger equation (NSE) 

idta(x, t) = -iw(2)a^^a(x, t) -f x|a(a;, i)pa(x, t), (1) 

(i, propagation variable; x, "transverse" coordinate; 
cj^^^ second-order dispersion or diffraction constant; x, 
nonlinear-coupling constant). Bright temporal solitons 
can be formed either in focusing media with anomalous 
dispersion (x < 0, w^^^ > 0) or in defocusing media with 
normal dispersion (x>0, a;'-^-'<0). Spatial solitons al- 
ways require a focusing nonlinearity. Instead of working 
in x-space, we can also turn to the w-space. 
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In what follows the coordinates in the x- and tj-domains 
are scaled by the initial soliton width xq and the spec- 
tral width Wo = l/s^o respectively. Propagation dis- 
tances are measured in dispersion (diffraction) lengths 
td = a;^/|w(2)|. The NSE @) can be solved by the inverse 
scattering method In particular, the bell-shaped 

fundamental soliton is extremely stable and can prop- 
agate over distances far beyond the limit given by the 
dispersion- (diffraction-) assisted spreading. Scaling the 
power of the fundamental soliton by a factor of A^^, A^ = 
2,3, . . ., yields so-called A^-bound solitons |2^]. The two 
classes of solutions of the NSE have been realized exper- 
imentally (see, e.g., |2^-|2|]). 

In quantum theory, the complex field amplitude a{x,t) 
and a* (a:, t) become field operator a{x, t) and a){x, t), re- 
spectively, with 



[a{x , t) , a' {x' , t)\ =S{x — x'), 



(3) 



and the NSE becomes an operator-valued equation 
[^ , pT| . In order to solve the quantum mechanical prob- 
lem, it is commonly assumed that the initial state is a 
multimodc coherent state such that (a) ~ A^sech(a;/a:o) 
is the initial shape of the classical A^-bound-soliton so- 
lution. For not too large propagation distances the state 
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then evolves into a multimode squeezed state of Gaus- 
sian type, which can be calculated numerically within 
the framework of appropriate discretization . 

The nonlinear soliton motion leads to internal quan- 
tum correlations, which have been studied both experi- 
mentally [H and theoretically |||,||-|||). From the re- 
sults the question is suggested of whether or not it is 
possible to produce two entangled light beams by appro- 
priate partitioning of quantum solitons. Let us consider 
bipartition as sketched in Fig. |l|, which may be realized 
experimentally using spectral filtering [ p7| , p8| or spatial 
filtering ]3l] ]. Since the soliton quantum state remains 
pure during the propagation, the entanglement E of the 
bipartite system is given by the von Neumann entropy of 
a subsystem [O i.e.. 



E = Si = -Tri(/5iln/5i) = ^2 



(4) 



where pi is the density operator of the zth subsystem, 
and Tr^ means the trace with regard to the ith subsys- 
tem (i = l,2). 



point of partition 



subsystem 1 
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entanglement ? 



FIG. 1. Bipartition of optical solitons. Applying appropri- 
ate filtering techniques (see, e.g., ||3l| , ^7|j3^ ) , it can be realized 
in the a;-domain by partition of the near field of spatial soli- 
tons and in the tj-domain by partition of the far field of spatial 
solitons or spectral partition of temporal solitons. 



In order to calculate the entropy of a subsystem pre- 
pared in a Gaussian state, let us consider the x-domain 
and assume that the subsystem extends over a region S 
(the calculation in the w-domain is analogous). We intro- 
duce new bosonic operators 



bk = dx [p.k{x) Aa{x) + Vk{x) lS.a\x)\ 



(5) 



[Aa(x) = a{x) — {a{x))] and impose on the pk{x) and 
Vk [x) the condition that 
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For notational convenience, we have omitted the argu- 
ment t. From Eqs. (^) and (^) it then follows that the 
parameters ak and the functions pk(x) and Vk(x) solve 
the equations 



(ct/c + 5) 5kk' 
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C{x,x')* B{x,x') 
B{x,x'y C{x,x') 



p.k{x) Vk{x) 



where 



B{x,x') = (Aa(a;)Aa(x')), 

C{x, x') = (Aat(x)Aa(a;')> + \5{x - x'). 



(7) 



(8) 
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Note that the functions p,k{x) and Vk{x) must satisfy the 
conditions 



dx [pk{x) pI,{x) - Vk{x)vl,{x)] =6kk', (10) 



dx [pk{x) i^k'{x) - i^k{x) p.k'{x)] =0, (11) 



because of the commutation relations [bk,b]^,] —Skk' and 

[bk,bk'] — 0. With regard to the new bosonic operators, 
the Gaussian state (of the subsystem) is obviously a mul- 
timode thermal state, with ak being the mean photon 
number of the fcth mode. Hence, the entanglement E is 
given by 



E = ^ Sth{<Jk), 



with 



Sthicr) = In 



(12) 



(13) 



being the entropy of a single-mode thermal state. Note 
that when the S-region comprises the whole soliton pre- 
pared in a pure Gaussian state, then all the ak vanish 
and thus E = 0, as expected. 

Figure ^ presents examples of the entanglement ob- 
tained by symmetrical bipartition of the fundamental 
soliton and the 2- and 3-bound solitons in both the x- 
domain [solid lines in Fig. §(6)] and the cj-domain [solid 
lines in Fig. ^(c)] as a function of the propagation dis- 
tance. The figure clearly shows that highly entangled 
parties can be realized by bipartition of optical solitons 
produced in Kerr-like media. The (numerical) calcula- 
tions have shown that the maximum entanglement is typ- 
ically observed for symmetrical bipartition (i.e., at a; = 
or Lu = 0). For the system considered in the figure, an 
asymmetrical bipartition of the 3-bound soliton in the x- 
domain at ^ 0.4x0 can improve the entanglement by 
5.3% for certain propagation distances. 
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The entanglement realized by bipartition of the fun- 
damental soliton increases monotonously with the prop- 
agation distance. The maxima and minima of the en- 
tanglement realized by bipartition in the x-domain of 2- 
and 3-bound solitons follow quite exactly the periodic 
change between soliton compression and expansion [com- 
pare Figs. ||(a) and ||(&)]. In the cj-domain this correla- 
tion is less pronounced [compare Figs. ||(a) and ||(c)]. It 
is worth noting that the amount of entanglement that 
can be achieved increases with the order parameter TV. 

Formally, the sum in Eq. ( |l2|) runs over all modes. In 
fact, only a few modes substantially contribute to the en- 
tanglement, as it is seen from Figures ^ andjH (the modes 
are numbered such that the single-mode contribution to 
the entanglement decreases with increasing mode index) . 
In particular, it is seen that the number of relevant modes 
is just given by the soliton-order parameter N. A possible 
reason for this fact can be seen in the multicomponent 
structure of the classical iV-soliton solution as discussed 
in @. 

In summary, the results offer novel possibilities of us- 
ing optical solitons as sources of entangled light beams. 
Although bipartition of solitons again yields multimode 
objects, only a few modes are involved in the entangle- 
ment. In particular, fundamental solitons are suited for 
generation of entangled light like 2-mode squeezed light, 
which is an alternative to the standard x^^^ sources. 
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FIG. 2. Dependence on the propagation distance of the entanglement realized by symmetrical bipartition of the fundamental 
soliton {N — 1) and Ai'-bound solitons (A'' = 2,3) in the a;-domain (b) and the ctJ-domain (c). For comparison, the soliton 
mid-intensity in the x-domain as a function of the propagation distance is shown (a). The contributions to the entanglement 
of the relevant modes are shown for the fundamental soliton and the 2-bound soliton. The marked propagation distances are: 
t4 — ^ td (soliton period), t2 = \t4 (distance of compression for 2-bound soliton), ti — jti and fa = 1 14 (distances of compression 
for 3-bound soliton). In the numerical calculations, which are performed on a grid of 256 points with a discretization step 
Aa; — 0.05 xo, it is assumed that the total photon numbers of the solitons are 2N^n, where n=10^ (for details, see [0,||]). 
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FIG. 3. The contributions of the relevant modes to the entanglement realized by symmetrical bipartition of the 3-bound 
soliton in the a;-domain (a) and the CLi-domain (6). The soliton is the same as in Fig. ^. 
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